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We show how the deposition of laser energy induces propulsion and strong deformation of an ab¬ 
sorbing liquid body. Combining high speed with stroboscopic imaging, we observe that a millimeter¬ 
sized dyed water drop hit by a millijoule nanosecond laser pulse propels forward at several meters 
per second and deforms until it eventually fragments. The drop motion results from the recoil 
momentum imparted at the drop surface by water vaporization. We measure the propulsion speed 
and the time-deformation law of the drop, complemented by boundary-integral simulations. The 
drop propulsion and shaping are explained in terms of the laser-pulse energy, the drop size, and the 
liquid properties. These findings are, for instance, crucial for the generation of extreme ultraviolet 
light in nanolithography machines. 


INTRODUCTION 

Laser-induced phase change in liquids can lead to a vi¬ 
olent response: deformation and disruption of the liquid 
body followed by the ejection of matter. The complete 
vaporization or even explosion of micrometer-sized drops 
can result from the linear absorption of laser energy [1- 
3]. Self-focussing and dielectric breakdown may lead to 
plasma formation in transparent drops [4-7]. Laser im¬ 
pact has been used to generate liquid motion by vaporiza¬ 
tion or plasma formation in confined geometries [8-10], 
sessile drops [11], and biological matter [12-14]. 

Here, we show how the absorption of laser energy by an 
unconfined liquid drop induces a rapid phase change (see 
Fig. 1), which in turn controls the propulsion, expansion, 
and fragmentation of the drop. A key application of the 
drop shaping by laser impact is found in laser-produced 
plasma light sources for extreme ultraviolet (EUV) nano¬ 
lithography [ 5, 16]. In these sources the shape, position, 
and stability of a liquid tin body directly affect the con¬ 
version efficiency of liquid tin to a plasma that emits 
EUV light. 

The detailed understanding of the hydrodynamic re¬ 
sponse of an opaque liquid drop to laser impact poses 
two fundamental challenges. First, one needs to resolve 
how momentum is transferred from the laser to the drop. 
Second, the subsequent deformation dynamics and frag¬ 
mentation of the drop after impact have to be quanti¬ 
fied. Although drop impact onto a solid substrate has 
been studied thoroughly (for a selection, see e.g. [17- 
23]), no consensus on the deformation dynamics has yet 
been reached and only few studies [18, 19, 24-26] focused 
on the fragmentation. 
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FIG. 1. (Color online) Laser pulses (A = 532 nm) impact¬ 
ing from the left on magenta-dyed water drops of radius 
R 0 — 0.9 mm. Images are taken 30 ps after impact with a 
color camera and diffusive backlight illumination, (a) White 
plasma glow [27] and violent ablation from the drop induced 
by a focused laser beam, (b) Fluorescence of the dye and ab¬ 
lation at the drop surface due to local boiling induced by a 
uniform laser irradiation. 


EXPERIMENTAL METHODS 

Our model system consists of a highly-absorbing drop 
that is hit by a pulsed laser beam. In Fig. 2 an overview of 
the experimental setup is shown. The drop detaches from 
a capillary, falls, and relaxes to a spherical shape with 
radius Ro = 0.9 mm. While it falls down, the drop masks 
a photodiode that then generates a reference trigger for 
the pulsed laser, cameras, and light source. The e^-axis 
of the laser beam is aligned orthogonally to the e y -axis 
defined by the falling drop and the e x -axis of the imaging 
optics. The xz-plane in which the laser beam propagates 
is below the plane of the trigger laser and the pinch-off 
point at the capillary tube. 

The drop consists of dyed water with a density p = 
998 kg/m 3 and surface tension 7 = 72mN/m assumed 
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to be equal to the properties of pure water. The typical 
penetration depth of the laser light into the dyed drop is 
S ~ 10 pm <C Ro [ 10 ], which ensures that the laser energy 
is absorbed in a thin layer close to the drop surface. The 
laser-pulse energy is varied between 0 and 120 mJ by an 
optical attenuator based on a half-wave plate and a polar¬ 
izing beam-splitter. The relation between the laser-pulse 
energy at the drop location and the settings of the atten¬ 
uator is determined in separate measurements, for which 
the top beam dump shown in Fig. 2 is replaced by an en¬ 
ergy meter. A focusing lens decreases the beam diameter 
to twice the drop size in order to achieve a uniform but 
high-intensity illumination of the drop. To ensure the 
drop is placed at the center of the laser beam, the drop 
position is optimized such that the drop-shape evolution 
is axi-symmetric with respect to e z and the propulsion 
speed is maximum. 

The energy E that is actually absorbed by the drop 
is computed from a beam-profile measurement and ray¬ 
tracing [28] . The typical beam fluence 1 J/cm 2 is well 
below the dielectric breakdown and self-focusing thresh¬ 
olds reported for water with focussed nanosecond laser 
pulses [ 8 , 27]. Consistently, we observe a plasma only 
when the laser beam is tightly focused inside the drop 
(Fig. la, see also [ 2, 13]), but not for a uniform irradia¬ 
tion (Fig. lb). To visualize the wavelengths in the visible 
spectrum that are emitted by the drop shown in Fig. 1 
we use a magenta-colored ink as a dye and a color cam¬ 
era. For all experiments leading to quantitative results 
we use a black-colored ink to suppress fluorescence. 

The post-impact dynamics of the drop (Fig. 3a) is ob¬ 
served from a side view (e^-axis in Fig. 2) with a long¬ 
distance microscope, a high-speed camera operated at a 
frame rate of 20 000 frames per second and a continuous 
light source. Detailed information in the first microsec¬ 
onds after impact is obtained by operating a camera in 
stroboscopic mode with a flash lamp that delivers a high- 
intensity light pulse of 8 ns (Fig. 3b). We record strobo¬ 
scopic videos by performing a single impact experiment 
per video frame while changing the time delay between 
the laser impact and the pulsed light source. For both 
cameras used the size of the field of view is 16x10 mm 2 , 
which yields a pixel resolution of 16 pm per pixel. 

RESULTS & INTERPRETATION 

The drop dynamics for different pulse energies is shown 
in Fig. 3. On impact, the surface of the drop hit by the 
laser emits a shock wave into the air (Fig. 3b). The shock 
wave is followed by the ejection of a mist cloud of small 
drops that is visible as a gray-to-black haze in the images 
and persists for several microseconds. Subsequently, the 
mist is expelled while the drop propels in the opposite di¬ 
rection (Fig. 3a). At the same time the drop flattens and 
expands in the radial direction before it either retracts, 



FIG. 2. (Color online) Sketch of the experimental setup in 
top view. The drop (black/magenta ink, IJC-5900/5920 by 
Sensient) is generated with a capillary tube connected to a 
syringe pump (PHD2000 by Harvard Apparatus, not shown). 
A continuous-wave laser (CWL, CPS196 by Thorlabs) and a 
photo diode (PD, PDA36A by Thorlabs) serve as light barrier 
to trigger on the falling drop. The lab equipment is synchro¬ 
nized by a high-precision pulse delay generator (BNC575 
by Berkeley Nucleonics) according to the indicated signal 
path. The pulsed laser is a frequency-doubled NdiYAG laser 
(Evergreen 140 by Quantel) with a pulse duration r p = 10 ns 
emitting at a wavelength A = 532 nm. Attenuation of the laser 
energy is accomplished by a zero-order half-wave plate (A/2), 
a polarizing beam splitter (PBS), and a beam dump (BD). 
The laser-pulse energy is measured by an energy meter (EM, 
QE12 by gentec-eo). The circle and arrow symbols along the 
optical axis respectively indicate the S- and P-component of 
the linearly-polarized laser beam that is focused by a plano¬ 
convex lens with a focal length of / = 125 mm. Side-view 
images (yz- plane) are taken with a long-distance microscope 
(LDM, 12x Zoom by Navitar), a high-speed camera (FAST- 
CAM SA-X2 by Photron), and a continuous light source 
(LS-M352A by Sumita). Stroboscopic images are acquired by 
a CCD camera (PCO1300 by PCO AG) combined with a 
light source (NANOLITE KL-K by HSPS) that delivers a 
high-intensity light pulse of 8 ns. A notch filter protects the 
imaging equipment from scattered laser light. 


for low pulse energy, or fragments, for large energy. 

We quantify the drop motion by measuring the dis¬ 
placement Z(t) of the drop center-of-mass and the drop 
radius R(t) (defined in Fig. 3a) for the first millisec¬ 
onds after impact. As Fig. 4a shows, the drop is pro¬ 
pelled at a constant speed U that increases with increas¬ 
ing pulse energy up to 2.0 m/s. The accompanying de¬ 
formation of the drop occurs on the inertial time-scale 
T[ = Ro/U ~ 10 -4 to 10 -3 s (Fig. 4b) and is eventually 
slowed down by surface tension on the capillary time- 
scale r c = yj pRq /7 = 3.5 ms. Both the initial deforma¬ 
tion rate r i _1 and the maximal extension R max increase 
with increasing pulse energy. We emphasize the clear 
separation of time-scales 

Tp < T e < 7i < T c 
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FIG. 3. (Color online) Side-view of a dyed water drop with initial radius Ro = 0.9 mm hit at t = 0 by a laser pulse propagating 
from left to right (e z ). (a) Drop shape dynamics for pulse energies increasing from bottom to top. E is the energy that is 
actually absorbed by the drop, We is the Weber number of the propelled drop (see text). The images are taken at a frame 
rate of 20.000 frames per second (r c =3.5ms). As the laser ablates the front of the drop a mist cloud is ejected backward (-e z ) 
while the remainder of the drop is propelled forward (e z ) and expands radially (e r ). For small E the drop retracts after the 
initial expansion and no break-up occurs. For E — 24 mJ the edge destabilizes before it retracts and the drop fragments, (b) 
Close-up view of the drop surface for E — 24 mJ revealing the shock wave in the air and the mist cloud development at early 
times (increasing from top to bottom). 


between the successive steps, namely, the laser pulse, the 
ejection of matter (on time scale r e ~ 10 -5 s), the initial 
deformation of the drop, and its capillary retraction. 

To explain the relation between the drop propulsion 
speed, the radial expansion, and the laser energy one 
needs to understand the mechanism that propels the 
drop. Surely, both the optical radiation pressure from the 
laser and the thermal radiation pressure caused by the 
heating of the drop surface are insignificant [29, 30]. The 
motion actually results from the recoil due to the par¬ 
tial vaporization of the drop: since the highly-absorbent 
dye ensures that the laser energy is absorbed in a super¬ 
ficial layer on one side of the drop, the vapor expulsion 
is mainly unidirectional and consequently transfers mo¬ 
mentum to the remainder of the drop. 

The light energy is absorbed by a liquid mass ~ pR^S 
set by the penetration depth of the laser. On the time 
scale T e , both diffusive and radiative heat transfers are 
negligible (the thermal diffusion length is much smaller 
than S [31]). Since the beam profile is flat, and neither 
the focusing due to the drop interface curvature nor non¬ 
linear optical effects (self-focusing or electric breakdown) 
are significant, we consider the energy deposition in the 
superficial layer to be close to uniform. This energy is 
sufficient to heat the liquid from the ambient tempera¬ 


ture T 0 = 293 K to the boiling temperature T b ~ 393 K, 
but not to vaporize all of it: only a certain fraction (3 
actually vaporizes. The energy balance therefore reads 
E ~ pR%5[c v (T h -T 0 )+(3AH], where c v = 4.0kJ/(kgK) 
and AH = 2.25 MJ/kg are, respectively, the specific heat 
capacity and latent heat of vaporization of the liquid. 

In all our experiments a mist cloud is observed, which 
is a clear signal of a local boiling of the drop. We there¬ 
fore assume that to get propulsion, a threshold energy 
T t h ~ pRqSc v (T b — Tq) ~ 3mJ has to be absorbed by 
the superficial layer to heat the liquid to the boiling point, 
which is in good agreement with the threshold for propul¬ 
sion observed in our experiments (Fig. 5a). Any addi¬ 
tional energy deposited in the superficial layer is used to 
vaporize a mass of liquid m ~ /SpR^S ~ (E — E^)/AH. 
An upper limit for the proposed scaling is given by 
T/T t h ~ 1 + Ai7/[c v (T b — To)] ~ 8, in which case the ab¬ 
sorbed energy is sufficient to evaporate the entire heated 
liquid layer (i.e. /3 = 1). Any increase in E beyond this 
point would lead to a superheated or even a critical phase, 
in which case the opaque mist cloud would not be ob¬ 
served [13]. 

For 0 < (3 < 1, which is the case of our experiments, 
the remaining part of the heated layer that is not va¬ 
porized is expelled as a mist of small drops. We assume 
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FIG. 4. (Color online) Center-of-mass displacement Z (along 
e z , see Fig. 3) (a) and radial expansion R (b) as a function of 
time for different absorbed energies. The corresponding image 
sequences are shown in Fig. 3. Each point is averaged over 
two experiments and the shaded area indicates the difference 
between the two. The apparent acceleration in Z for t/r c < 
0.2 is an artifact of the method used to determine the center- 
of-mass position [28]. The large deviation in R for E = 24 mJ 
illustrates the statistical nature of the fragmentation. For 
E — 15 mJ events of drop ejection from the edge are visible 
at t/r c — 0.24, 0.36 and 0.63. 



b) 



FIG. 5. (Color online) (a) Propulsion speed of the drop as a 
function of the absorbed laser energy. Each point represents 
at least four experiments, the error bars indicate the standard 
deviation. The solid line is equation (2) with a prefactor of 
0.4. (b) Maximal relative expansion Rma^/Ro — 1 for individ¬ 
ual experiments as a function of the Weber number in linear 
and logarithmic (inset) scales. The solid line is equation (3) 
with a prefactor of 0.6. For large Weber numbers a saturation 
is observed due to the fragmentation of the sheet. 


that the liquid vaporizes at Xb and that the vapor is ex¬ 
pelled at the thermal speed u = ^/k^T^/p ~ 400 m/s, 
where k & ~ 1.38 x 10 -23 J/K is the Boltzmann constant 
and fj, = 2.99 x 10“ 26 kg is the molecular mass of water. 
This expelled vapor propels the remainder of the drop. 
Momentum conservation mu = pR^U yields 


U 


E — E th 
pR%ah u ’ 


( 2 ) 


that is, an increase in U proportional to that in E. Figure 
5a shows that this scaling argument, with a prefactor of 
0.4, is in good agreement with our experimental data. 

With a description of the propulsion at hand, we now 
turn to the drop deformation. The expansion dynamics is 
directly affected by surface tension, which promotes the 
retraction and possibly the fragmentation of the drop. 
The key parameter describing the expansion is therefore 
the Weber number of the motion induced by the laser 
We = pRoU 2 /y, which compares the drop displacement 
kinetic energy to its surface energy. In our experiments 
1 < We < 60. The impulsive acceleration of our drop 
from 0 to U is similar to the impulsive stop of a drop 
impacting a solid with velocity U. We therefore use the 
momentum-based scaling derived by [17, 19] for drop im¬ 


pact on solids to express the maximal radial expansion 


^max Ro 

Ro 


We 1/2 


I pR 0 u 2 E - E th 

7 m Air 


(3) 


in which the expression in terms of E directly comes from 
(2). Expression (3), with a prefactor of 0.6, is in good 
agreement with our experimental data up to We ^ 40, 
when the drop starts fragmenting and the maximum ex¬ 
pansion saturates (see Fig. 5b). The scaling (3) has al¬ 
ready been observed for drop impact onto solid substrates 
with negligible friction [19]. The present setup is however 
fundamentally different since, as mentioned above, the 
typical impact timescale, during which the drop acceler¬ 
ates, is decoupled from the inertial timescale: r e <C R. 


NUMERICAL RESULTS 

To confirm that the interaction of the laser pulse with 
the drop can be modeled as a short recoil-pressure pulse 
exerted on the drop surface, we perform boundary inte¬ 
gral (BI) simulations [32-34]. We assume that the flow 
inside the drop is inviscid, irrotational, and incompress¬ 
ible, and solve the resulting Laplace equation for the flow 
potential. The method assumes axi-symmetry and there- 
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fore cannot be used to study the eventual fragmentation 
of the drop, but it does capture the initial phase of the 
drop deformation. 

The laser pulse is modeled by applying a pressure 
boundary condition at the drop surface for a time du¬ 
ration r e < Ti. We use a Gaussian pressure profile with 
a length-scale based on the measured laser-beam profile 
and a pressure scale set to match the propulsion velocity 
observed in the experiment. This pressure scale is pre¬ 
scribed by the momentum conservation pR^r e r\j pRlU 
(the prefactor can be obtained analytically [35]). From 

(2) this recoil pressure can readily be expressed in terms 
of the absorbed energy. 

The numerical drop shape evolution is shown in Fig. 6. 
It illustrates the added value of the simulations: not only 
the two-dimensional projection of the drop shape, but 
also the spatial and temporal evolution of the sheet thick¬ 
ness can be extracted, which is crucial when it comes to 
study fragmentation [L9]. Moreover, Fig. 6 shows that 
the BI model quantitatively predicts the radial drop ex¬ 
pansion observed for different Weber numbers. This con¬ 
firms that a pressure pulse applied at the drop surface for 
a time much shorter than the hydrodynamic time scales 
(Ti and r c ) is indeed sufficient to describe the hydrody¬ 
namic response of a drop to the impact of a laser pulse. 


CONCLUSIONS 

We have shown that an opaque free-falling drop hit 
by a laser pulse propels and expands until fragmentation 
occurs. In the present case the laser energy is absorbed 
in a superficial layer of the drop such that the deposited 
energy per unit mass E/pR^S ~ 0.1 to lMJ/kg is com¬ 
parable to the specific latent heat of vaporization. As a 
consequence, drop motion is induced by the recoil due 
to vaporization on the face of the drop that is hit by 
the laser. This results in a propulsion speed and a max¬ 
imal radius of expansion that are both proportional to 
the pulse energy. The expansion dynamics is limited by 
surface tension and is similar to that of a drop impacting 
a solid, although with a laser pulse momentum transfer 
takes place on a much shorter time scale. Laser-induced 
drop fragmentation and the influence of the beam fo¬ 
cussing require detailed studies and are left for future 
work [35]. All results reported here should transpose di¬ 
rectly to the shaping of liquid tin drops in EUV light 
sources. In a regime where a plasma is generated the 
propulsion mechanism may change, however, the Weber 
number remains the key parameter governing the hydro- 
dynamic response. 

We thank Chris Lee, Andrea Prosperetti, and Guil¬ 
laume Lajoinie for fruitful discussions. This work is part 
of an Industrial Partnership Programme of the Founda¬ 
tion for Fundamental Research on Matter (FOM), which 



FIG. 6. (Color online) Radial expansion of the drop: exper¬ 
iments (markers) and BI simulations (solid lines). The cor¬ 
responding image sequences are shown in Fig. 3. The inset 
shows the drop shape evolution from BI for We = 60 with 
an exaggerated center-of-mass displacement to separate the 
successive shapes. The simulations are stopped when the lo¬ 
cal sheet thickness becomes too thin to resolve the dynamics 
accurately. 
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